FALL3D: A computational model for transport and deposition of volcanic ash

FALL3D is a 3-D time-dependent Eulerian model for the transport and deposition of volcanic ashes and lapilli. The model solves the advection–diffusion–sedimentation (ADS) equation on a structured terrain-following grid using a second-order finite differences (FD) explicit scheme. Different parameterizations for the eddy diffusivity tensor and for the particle terminal settling velocities can be used. The code, written in FORTRAN 90, is available in both serial and parallel versions for Windows and Unix/Linux/Mac X operating systems (OS). A series of preand post-process utility programs and OS-dependent scripts to launch them are also included in the FALL3D distribution package. Although the model has been designed to forecast volcanic ash concentration in the atmosphere and ash loading at ground, it can also be used to model the transport of any kind of airborne solid particles. The model inputs are meteorological data, topography, grain-size distribution, shape and density of particles, and mass rate of particle injected into the atmosphere. Optionally, FALL3D can be coupled with the output of the meteorological processor CALMET, a diagnostic model which generates 3-D time-dependent zero-divergence wind fields from mesoscale forecasts incorporating local terrain effects. The FALL3D model can be a tool for short-term ash deposition forecasting and for volcanic fallout hazard assessment. As an example, an application to the 22 July 1998 Etna eruption is also presented. & 2008 Elsevier Ltd. All rights reserved.